// Ninf-G-XX, ver. 0.1
// Copyright 2011 National Institute of Advanced Industrial Science and Technology
// 
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
// 
//     http://www.apache.org/licenses/LICENSE-2.0
// 
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.
#include <stdio.h>

void matmul(int n, double * A, double * B, double * C){
  int i, j, k;
  printf ("n = %d\n", n);
  for (i = 0; i < n; i++){
    for (j = 0; j < n; j++){
      double sum = 0;
      for (k = 0; k < n; k++)
	sum +=  *(A + i + k * n) * *(B + k + j * n);
      *(C + i + j * n) = sum;
    }
  }
}

void matsub(int n, double * A, double * B, double * C){
  int i, j, k;
  printf ("n = %d\n", n);
  for (i = 0; i < n; i++){
    for (j = 0; j < n; j++){
      *(C + i + j * n) = *(A + i + j * n) - *(B + i + j * n);
    }
  }
}
